Parameterized optimized effective potential for atoms 
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The optimized effective potential equations for atoms have been solved by parameterizing the 
potential. The expansion is tailored to fulfill the known asymptotic behavior of the effective potential 
at both short and long distances. Both single configuration and multi configuration trial wave 
functions are implemented. Applications to several atomic systems are presented improving previous 
works. The results here obtained are very close to those calculated in either the Hartree-Fock and 
the multi configurational Hartree-Fock framework. 
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I. INTRODUCTION 

The Hartree-Fock (HF) method is the best approxi- 
mation to the atomic and molecular problem in the inde- 
pendent particle model with a single configuration. Orig- 
inally it was formulated in terms of a single Slater de- 
terminant but in current applications a single configura- 
tion, with one or more Slater determinants coupled in 
the LS scheme, is used, see e.g. P]. An alternative to 
the HF equations, always within the independent par- 
ticle model, is given by the Optimized Effective Poten- 
tial (OEP) method proposed by Sharp and Horton I3| as 
a variational alternative to the Slater's simplified treat- 
ment of the exchange term in the HF equations based on 
the averaging in the occupied orbitals [Sj . 

In the OEP method an additional constraint is im- 
posed into the variational problem: the orbitals must 
satisfy a single-particle Schrodinger equation with a cer- 
tain local potential, the same for all the electrons (with 
the same spin). The expectation value of the hamilto- 
nian of the N-electron system becomes a functional of 
such a local potential. The effective potential is then 
varied to minimize the total energy. This gives rise to a 
linear integral equation in the effective potential 0, H 
whose solution gives the optimized effective potential. 
The wave function of the system is then constructed from 
the single-particle wave functions which are eigenfunc- 
tions of the optimized effective potential. 

The OEP method has been used in connection with the 
Kohn-Sham density functional theory 0, S S • Accu- 
rate spin-polarized exchange-only Kohn-Shan potential 
has been constructed with the OEP method |j9|. The lo- 
cal exchange potential obtained from the OEP method 
has many of the analytical features of the exact Kohn- 
Sham potential, and it has been recognized as the exact 
implementation of exchange-only density functional the- 
ory |Tp| . For example the HF potential can support very 
few unoccupied excited states because of its exponential 
fall off and therefore it is not a good starting point to 
describe excited states. On the other hand as the OEP 
potentials presents the proper long range behavior it pro- 
vides better excitation energies. 

The use and application of the original form of the 



OEP was hindered because of the computational diffi- 
culties posed by the OEP integral equation. Initially, 
a numerical grid method was employed by Talman and 
Shadwick P| to solve this equation and further more ac- 
curate calculations for atoms have been carried out by 
using a more refined mesh HI [13. Also some approxi- 
mations to the OEP equations were proposed ISj. Very 
recently an iterative scheme based on solving the par- 
tial differential equations satisfied by the orbitals shift 
that satisfy exactly the Krieger, Li and lafrate approxi- 
mation has been devised An alternative methodol- 
ogy to board this problem uses a finite basis set expan- 
sion and /or a naly tical forms for the effective potential 
lilllllallllialllllSIIlIll- The parameterization 
of the effective potential simplifies greatly the numerical 
problem of solving the integral equations 0, 0, . 
The explicit form or the parameterization is proposed by 
using different arguments ^16, .19t , .23i] . The parameters 
can been fixed either by using the variational method or 
by matching the one-body eigenvalues of the correspond- 
ing spin-averaged Dirac equation to the experimentally 
observed one-electron ionization energies \u\ . 

The aim of this work is to obtain the energy and other 
properties for atoms in the parameterized OEP approxi- 
mation. The analytical form used for the optimized po- 
tential has been used previously 0, ^3 but with sim- 
pler expansions. In this work we increase the number 
of free parameters in the effective potential until conver- 
gence is reached. In this way we improve previous results 
obtained not only within the parameterized OEP scheme 
but also by solving numerically the integral equation. We 
have studied the ground state of the atoms Li to Ar in the 
LS coupling. Finally the OEP approximation is extended 
here to multi configuration wave functions in the same 
spirit of the Multi Configuration Hartree Fock method. 
The binding energy, single particle energies and exchange 
energies are obtained and compared systematically with 
their Hartree-Fock counterparts. The best results for the 
energy of the present work are roughly 1 mhartree above 
the best self consistent field results showing the good per- 
formance of the local potential approximation. 

The structure of this work is as follows. In Section UTI 
we show in detail the parameterization and the algorithm 
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used. The results obtained are reported and discussed 
in Section 11111 The conclusions and perspectives of this 
work can be found in Section llVI Atomic units are used 
throughout. 



II. PARAMETERIZED OPTIMIZED 
EFFECTIVE POTENTIAL 

The single-particle wave functions used to construct 
the Slater determinants of a given single configuration 
are the eigenf unctions of the so called effective potential, 
that in this work is taken to be central, 



+ KM j ^xAr) = eA0A,.(r) (1) 



where A and tr stand for the spatial and spin quantum 
numbers, respectively. 

The total energy of the system is therefore a functional 
of this single particle potential. The minimum condition 
on this potential leads to the OEP equations |3 fixing 
the best effective potential. The energy so obtained will 
be an upper bound to the exact one and it will be above 
the HF value. It is also worth pointing out that the role 
of the effective potential here is just an auxiliary function 
used to calculate the orbitals in the wave function. 

The parameterization chosen in this work for the effec- 
tive potential is 

/ S "p \ 

Ve{r) = ^- [Z - N +1 + {N -l)Y,Y. CKv^^e-^--^^ 

^ \ p=Ofe=l / 

(2) 
(3) 



with the condition 



fc=i 

imposed in order to match the correct short-range be- 
havior of the potential. This constraint makes that we 
must deal with, at least, two terms with p = 0, i.e. 
no > 2. With respect to the long range asymptotic be- 
havior, this functional form is such that the potential 
goes as — (Z — iV -f l)/r for large electron-nucleus dis- 
tances. The number of basis functions used to expand 
the effective potential is incremented systematically until 
convergence is reached. As we shall see later the impor- 
tance of considering higher values of p increases with the 
number of electrons. Finally it is worth mentioning here 
that regular behavior close to the nucleus has shown to 
be relevant in local effective potential theories p^ . 
The algorithm is the following 

1. Pick the initial values of the effective potential pa- 
rameters, Cfc.p and /3fc.p such that the short and long 
range asymptotic behaviors are satisfied. 

2. Solve the one-electron Schrodinger equation for 
each occupied orbital. 



3. Use these orbitals to build up the N-electron wave 
function and calculate the expectation value of the 
Hamiltonian. 

4. Optimize the total energy with respect to the free 
parameters (c^^p, /3fc,p) with the constraint given by 
equation 13 

5. Increase the size of the expansion of the effective 
potential until convergence. 

The single-particle Schrodinger equation ^ in step 3 
is solved by expanding the radial orbitals in a basis set 
of Slater type functions as in the Roothaan-Hartree-Fock 
method. In particular we have used the same size of the 
basis set as in reference [S^: the s-type single particle 
orbitals are expanded by using 6 basis functions for the 
Li to Ne atoms and 8 basis functions for the Na to Ar 
atoms. The radial part of p-type single particle orbitals 
is developed as the sum of 4 basis functions for the B 
to Ne atoms, 5 basis functions for Na and Mg and 8 ba- 
sis functions for the atoms Al to Ar. By doing this all 
the different matrix elements can be computed analyti- 
cally. The total energy has been minimized with respect 
to the free parameters of the effective potential by using 
a SIMPLEX 26] algorithm. 

Once obtained the single particle wave functions we 
can calculate some other quantities of interest such as 
the exchange energy 



r — r 



(4) 

and the Hartree-Fock single-particle energies, , ob- 
tained from the expectation value 



(5) 



where I, J and /C are the usual single particle, direct 
and exchange terms calculated starting from the eigen- 
functions of the effective potential. The single particle 
energies e^^ of equation does not coincide with the 
eigenvalues t\ of equation ^ except for the highest oc- 
cupied orbital. This was proven in reference |27| within 
a framework of spin-polarized and arbitrary exchange- 
correlation functionals. In the scheme of the present 
study this condition states that = , where Xh 
stands for the highest occupied level. The fulfillment of 
this condition has been used previously, e.g. P, T?] to 
asses the accuracy in the solution of the OEP equations. 

In general the optimization procedure is very stable for 
both the total energy and the Hartree-Fock values e^^. 
However this is not the case for the eigenvalue e\, as has 
been noted previously 0, |3| by using a numerical grid 
method. For this reason, and because the exact solution 
of the OEP must also satisfy the virial relation and the 
exchange-only virial relation |28l l29l | we have minimized 
the quantity 



{H)+ I CA. 



— E" 



(6) 
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where EJ^^ is given by 

E-"/ ^- j ATp{r)r-VV^{T) (7) 

and is the exchange potential. By doing this no signif- 
icant changes are found for both the total energy and the 
Hartree-Fock eigenvalues. One should expect that this 
procedure gives rise to a better description of the asymp- 
totic region of the optimized effective potential. As it has 
been previously pointed out these two conditions are 
in some sense complementary. The main contributions to 
the quantities involved in the exchange only virial condi- 
tion arise from the internal region of the atoms whereas 
the highest energy eigenvalue is governed by the outer 
region. By including these two conditions in the energy 
functional we observed a better and more stable conver- 
gence in the free parameters as compared to an uncon- 
strained minimization of {H). With the basis set used 
here both conditions are satisfied within one part in 10~^ 
hartree. 

The method can be generalized straightforwardly to 
deal with a Multi Configuration expansion. The starting 
trial wave function is written as a linear combination of m 
single configuration wave functions with the total orbital 
angular momentum and spin of the state under study. 
The hamiltonian is diagonalized in this set. The orbitals 
required to build the different Slater determinants are ob- 
tained as the eigenfunctions of the single particle effective 
potential containing several free parameters. The total 
energy is minimized with respect to those parameters as 
before. This more general trial wave function will provide 
not only a better description of the lowest energy state of 
a given symmetry but also a variational approximation to 
the excited states because of the Hylleraas-Undheim the- 
orem which states that the eigenvalues constitute upper 
bounds to the first m bound states. These type of wave 
functions have been recently used along with a correla- 
tion factor of Jastrow type to study some excited states of 
the berylhum atom of and its isoelectronic series (30i..31.j . 

III. RESULTS 

In Table m the results obtained by using different pa- 
rameterizations of the effective potential are analyzed. 
We show the values for the Ne atom which is representa- 
tive for the systems studied here. We compare with the 
numerical optimized effective potential results of Refs. 

and '33| and with the finite-basis set expansion of 
the effective potential of Refs. .19.. .22] . The approximate 
values of reference are also reported. The HF results 
are taken from |33l | and are the benchmark values for the 
total energy obtained from the OEP method. We also 
report the exchange energy, E^, and the eigenvalues e„; 
and ■ The notation used for the parameterization is 
(0"° , , 2"^ , . . .) where Up is the number of functions of 
the type e~'^'='p''r^'/r used in the expansion of the effective 
potentially 



The best results are obtained with the (0^, 1^, 2^) pa- 
rameterization, and this is the one used for the rest of 
the atoms. It can be seen that the value of the total en- 
ergy does not depend substantially on the basis set, but 
this is not the case of the other quantities. The use of 
bigger basis set sizes for the effective potential does not 
improve noticeably the energy for the atoms Li to Ar and 
it increases the computational effort. However for heav- 
ier atoms we have numerically checked that the rate of 
convergence can be substantially improved by including 
higher powers of r in the parameterization of the effective 
potential. 

The best results of this work improve previous ones 
obtained within the optimized effective potential scheme 
and it is only 0.57 mhartree above the Hartree-Fock 
result. The expectation values e^^ obtained with the 
OEP orbitals are in a very good agreement with the cor- 
responding single particle energies obtained within the 
Hartree-Fock framework. This is because of the fact, 
pointed out previously j^, that the single particle wave 
functions calculated from the OEP method are a good 
approximation to the Hartree-Fock orbitals. 

The ground state and exchange energies of the atoms 
Li to Ar as well as the single-particle expectation val- 
ues, e^f , and the single-particle OEP eigenvalues, e„i, 
are reported in tables and IIIII All of these quantities 
are compared with those obtained within the HF frame- 
work and with the numerical optimized effective poten- 
tial results of reference [s^- Within the Hartree-Fock 
framework the two quantities e^f and tni are the same 
so that we do not list them separately. For all the atoms 
the (0^,1^,2^) parameterization for the effective poten- 
tial has been used. 

As can be noted there is an appreciable improvement 
of the results, especially for the lightest atoms, which can 
be more easily seen in Figure^where we plot the relative 
error (in %) of our ground state energy with respect to 
the HF one as compared to the relative error for the nu- 
merical solution. In the present work the relative error is 
nearly constant for all the atoms considered. In principle 
one should expect a better energy coming from the nu- 
merical solution than from the parametrized one as the 
used in this work. 

Previous works have also reported energies calculated 
from a parametrized solution that are below the numer- 
ical solution 0, |23|. The reason for the better per- 
formance of the parametrized solution may be due to 
numerical inaccuracies in the numerical solution of the 
OEP equations due to the rather involved procedure of 
its solution. In an attempt to elucidate this fact we have 
carried out the following calculation. Starting from the 
tabulation of Aashamar et al. |32l | we have built a pa- 
rameterized potential of (0^,1^,2^) type by fitting the 
numerical values to that functional form. Then the fit- 
ted potential has been used in our code to determine the 
best orbitals within the effective potential approach of 
this work without changing the fitted potential. Here 
we show the results for the Ne atom, which are repre- 
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TABLE I: Convergence in the parameterization of the effective potential for the Ne atom. The results are compared with the 
Hartree-Fock values (HF), two different set of results obtained by solving numerically the optimized effective potential (NOEP) 
and with the Krieger, Li and lafrate approximation. (0"°, f"i,2"^, . . .) stands for Up functions of the type e~^'''P'^r^ /r in the 
expansion. 
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Figure 1, Sarsa et al. 

FIG. 1: Relative difference (in percent) between the OEP 
ground state energy of this work and the HF one for the atoms 
Li through Ar as compared with the same quantity calculated 
from the numerical results (NOEP) of ,32,1. 



effective manner. 

The difference between the numerical OEP energy and 
the resuh obtained here (—128.5465) is more accused. 
This discrepancy is due to the differences in the effective 
potential. In Figure [21 we plot the effective potential for 
the Ne atom and the differences with the numerical re- 
sults of ■ The results for Ne are representative for 
the rest of the atoms considered here. Notice that the 
differences are multiplied by 10 to better see them in the 
scale of the figure. 




sentative of the rest of the atoms analyzed. The total 
energy obtained in this way is —128.5457 hartree to be 
compared with —128.5455 hartree reported in Ref. [s^ . 
The difference between these two values is due to a better 
performance of the parameterized potential which allows 
one to work analytically. In particular the asymptotic 
behavior of the potential is taken into account in a more 



Figure 2, Sarsa et al. 

FIG. 2: Effective potential for Ne and the difference (multi- 
plied by 10) with the numerical results of ,32^1. 
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TABLE II: Energy and single particle energies for He to Ne atoms (OEP). The results are compared with the Hartree-Fock 
values (HF) obtained from reference |^ and with the numerical solution of the optimized effective potential equation of |33 |. 
the exchange energies for the closed shell atoms have been taken from 





Li('S) 


Be(iS) 


B(2p) 


C(^P) 


N(*S) 


0(»P) 


fCp) 


Ne(iS) 


E 
HF 
NOEP 
OEP 


-7.43273 

- i .4:0^4: 

-7.43261 


-14.57302 
-14.57291 


-24.52906 

94 ^97S 
-24.52874 


-37.68862 
-37.68774 


-54.40093 
-54.40029 


-74.80940 

74 8n7'i 

-74.80888 


-99.40935 
-99.40875 


-128.54710 

1 9S '\d'\'^ 

-128.54652 


Ex 
HF 
NOEP 
OEP 


-1.78119 
-1.78119 


-2.66692 
-2.666 
-2.66627 


-3.74759 
-3.74661 


-5.04930 
-5.05059 


-6.59707 
-6.59606 


-8.18189 
-8.18039 


-10.01105 
-10.01267 


-12.10835 
-12.105 
-12.10591 


Ml.' 

HF 
OEP 


-2.47774 
-2.47783 


-4.73267 
-4.73366 


-7.69534 
-7.69628 


-11.32552 
-11.32370 


-15.62906 
-15.63409 


-20.66866 
-20.66947 


-26.38276 
-26.37617 


-32.77244 
-32.77532 


eis 
NOEP 
OEP 


-2.08186 
-2.09495 


-4.12785 
-4.06912 


-6.91330 
-6.92455 


-10.35324 
-10.42413 


-14.46742 
-14.61828 


-19.21792 
-19.30917 


-24.66731 
-24.78314 


-30.81549 
-30.87242 


HF 

HF 
OEP 


-0.19632 
-0.19630 


-0.30927 
-0.30962 


-0.49471 
-0.49556 


-0.70563 
-0.70415 


-0.94532 
-0.94748 


-1.24432 
-1.24639 


-1.57254 
-1.56989 


-1.93039 
-1.93102 


NOEP 
OEP 


-0.19644 
-0.19630 


-0.30838 
-0.30962 


-0.52804 
-0.52968 


-0.74958 
-0.75608 


-0.99752 
-1.00990 


-1.19129 
-1.20004 


-1.43073 
-1.44285 


-1.71200 
-1.72224 


HF 

HF 
OEP 






-0.30986 
-0.30844 


-0.43334 
-0.43313 


-0.56759 
-0.56816 


-0.63191 
-0.63174 


-0.73002 
-0.72885 


-0.85041 
-0.85130 


e2p 
NOEP 
OEP 






-0.30976 
-0.30844 


-0.43067 
-0.43313 


-0.56322 
-0.56816 


-0.62933 
-0.63174 


-0.72479 
-0.72885 


-0.84571 
-0.85130 



In Figure El we plot, the function Vp{r) obtained form 
the effective potential as (see Eq. 

(8) 

We compare the best optimized effective potential ob- 
tained in this work with the numerical result calculated 
form the tabulated values of Ref. 32] , the fitted potential 
is included for the shake of completeness. 

Although the effective potential in the numerical ap- 
proach and the parameterized solution of this work is 
very similar, see Figure ^ the differences for Vp{r) be- 
come apparent. The main discrepancies are for r ^ 1.5 
au, with a faster decay to zero of the parameterized so- 
lution. The magnitude of this part of the effective po- 
tential is smaller than the total potential. Therefore, 
the discrepancies between the parametrized and numer- 
ical solutions are not very accused, in this case ~ 0.001 
mhartree. 

It is worth mentioning here that for Be, Ne, Mg and 
Ar atoms more accurate numerical results for highest 
occupied single particle eigenvalues have been reported, 
[E but the total energy is very similar to that of 
reference |3^ . 

To illustrate the performance of the Optimized Effec- 
tive Potential with a multiconfigurational trial wave func- 



tion we have applied it to a simple case. In Table HVI we 
study the ground and the first excited state of the 
beryllium atom using a two configuration wave function 
in both cases. We compare our results with the corre- 
sponding Multi Configuration Hartree Fock ones P, 01 
for the same states and configurations. For the ground 
state, where the 2s-2p near degeneracy plays an impor- 
tant role, we have expanded the wave function in terms of 
the configurations ls^2s^ and ls^2p^. We have also stud- 
ied the ^P state arising from the ls^2s2p configuration 
by considering the mixing with the ls^2p3d one. 

Thus for the ground state, the trial wave function is 
written as 

I*) = Ci|ls22s2;l S) + C2|ls22p2;l S) (9) 

and for the ls^2s2p ^P state the wave function is 

I*) = di\ls^2s2p;^ P) + d2\ls'^2p3d;^ P) (10) 

In Table HVl we show the values for the energy and for 
the coefficients Ck and dk obtain ed with the OEP method 
as compared with the corresponding MCHF ones [s^l ■ It 
is apparent the good agreement between the two sets 
of results for both states that illustrates that the OEP 
method not only provides a good value for the energy 
but also an adequate weight for any of the configurations 
involved in the corresponding state. For this case the 
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TABLE III: The same as in Table|n|for the atoms Na to Ar. 





Na(2S) 


Mg(iS) 


Al(2p) 


Si(3p) 


P(*S) 


S(3p) 


Cl(2p) 


Ar(iS) 


E 
HF 

NOEP 
OEP 


-161.85891 
-161.8565 
-161.85770 


-199.61464 
-199.6115 
-199.61303 


-241.87671 
-241.873 
-241.87415 


-288.85436 
-288.850 
-288.85160 


-340.71878 
-340.714 
-340.71537 


-397.50490 
-397.500 
-397.50155 


-459.48207 
-459.477 
-459.47854 


-526.81751 
-526.810 
-526.81405 


Ex 
HF 
NOEP 
(JEP 


-14.01752 
-14.01160 


-15.99429 
-15.988 
-15.98762 


-18.07225 
-18.06349 


-20.28350 
-20.27359 


-22.64091 
-22.62640 


-25.00614 

A onon 

-24.99093 


-27.51653 

OT /inooi^ 

-27.49826 


-30.18494 
-30.175 

on 1 'rnof^ 

-30.17936 


HF 

HF 
OEP 


-40.47850 
-40.47877 


-49.03174 
-49.03776 


-58.50103 
-58.50569 


-68.81246 
-68.81693 


-79.96971 
-79.97485 


-92.00445 
-92.01615 


-104.88442 
-104.89542 


-118.61035 
-118.60723 


eis 
NOEP 
OEP 


-38.04248 
-37.91728 


-46.32444 
-46.16018 


-55.5604 
-55.37384 


-65.6547 
-65.48096 


-76.59185 
-76.39372 


-88.34890 
-88.09753 


-100.9622 

1 r\r\ r rxr c\<-7 

-100.59527 


-144.4452 

1 1 .1 AT" '—rc\ r\ 

-114.45720 


HI-' 

HF 
OEP 


-2.79703 
-2.79384 


-3.76772 
-3.77025 


-4.91067 
-4.91136 


-6.15654 
-6.15701 


-7.51110 
-7.51443 


-9.00429 
-9.01030 


-10.60748 
-10.61285 


-12.32215 
-12.31768 


£23 

NOEP 

UEP 


-2.48096 
-2.22265 


-3.09619 
-3.11058 


-4.15912 
-4.17483 


-5.34036 
-5.355205 


-6.63449 
-6.64244 


-8.01012 
-7.97517 


-9.51097 
-9.54963 


-11.14504 
-11.19559 


HI.' 
£2p 

HF 
OEP 


-1.51814 
-1.51617 


-2.28223 
-2.28607 


-3.21830 
-3.22178 


-4.25605 
-4.25935 


-5.40096 
-5.40745 


-6.68251 
-6.69193 


-8.07223 
-8.07954 


-9.57146 
-9.56872 


NOEP 
OEP 


-1.19797 
-1.16901 


-1.85896 
-1.86617 


-2.73146 


-3.71956 
-3.73189 


-4.81814 

A 0000 r" 

-4.82825 


-5.99544 
-5.96260 


-7.29541 
-7.33117 


-8.72568 
-8.76408 


w 

£33 

HF 
OEP 


-0.18210 
-0.18172 


-0.25305 
-0.25359 


-0.39342 
-0.38663 


-0.53984 
-0.54090 


-0.69642 
-0.69896 


-0.87953 
-0.88305 


-1.07291 
-1.07617 


-1.27735 
-1.27597 


£33 
NOEP 
OEP 


-0.18211 
-0.18172 


-0.25206 
-0 25359 

\J ^ ^^^^ kJ 


-0.38859 
-0 38388 


-0.53937 
-0 53904 


-0.69510 
-0 69884 


-0.80812 
-0 80924 


-0.93838 
-0.94641 


-1.09280 
-1 09886 


£3p 

HF 
OEP 






-0.20995 
-0.21062 


-0.29711 
-0.29721 


-0.39171 
-0.39325 


-0.43737 
-0.43950 


-0.50640 
-0.50812 


-0.59102 
-0.58979 


esp 
NOEP 
OEP 






-0.20910 
-0.21062 


-0.29630 
-0.29721 


-0.38737 
-0.39325 


-0.43662 
-0.43950 


-0.50027 
-0.50812 


-0.58532 
-0.58979 



optimized effective potential method with muhi configu- 
ration trial wave function compares with the MCHF at 
the same level as the OEP method with the Hartree-Fock. 
More complex wave functions using a larger number of 
configurations have been used to study these and some 
other excited states of this atom and its isoelectronic se- 
ries [sol l3lj|. In those works, not fully optimized multi 
configuration wave function were used to build up more 
accurate explicitly correlated wave functions of Jastrow 
type. 



TABLE IV: Energy and coefficients for some multi configu- 
rational wave functions for the ground and the first excited 
state of ^P type for the beryllium atom. 



Is 


E 


Cl 


C2 


MCHF 


-14.61685 


0.95003 


0.31214 


OEP 


-14.61637 


0.95008 


0.31202 


^P 


E 




rf2 


MCHF 


-14.41156 


0.97524 


-0.22116 


OEP 


-14.41131 


0.97489 


-0.22270 



IV. CONCLUSIONS 

The optimized effective potential with parameterized 
potential has been used to study the ground state of 
the atoms from Li to Ar. Parameterized orbitals have 
been used to solve the corresponding single particle 



Schrodinger equation. The virial relation of Ghosh and 
Parr and Levy and Perdew involving the exchange energy 
and the exchange potential and a condition for the high- 
est energy occupied orbital of Krieger, Li and lafrate have 
been imposed. These are two analytically known condi- 
tions that the exact solution of the OEP equations must 
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Figure 3, Sarsa et al. 

FIG. 3: Vp(r) defined in Eq. Q obtained from tiie best 
parameterized effective potentiai of tfiis work, optimized, and 
the numericai resuits of Ref. numerical. Tfie anaiytical 

fit done to tiie numerical values is also plotted, fitted. 



fulfill. We have included them in our functional and a 



constrained search of the optimum parameters is carried 
out. As a result for the minimum both of them hold 
within 10~^ hartree and the minimum energy is not sub- 
stantially different to that obtained in an unconstrained 
minimization. The tail of the effective potential is ex- 
pected to be better reproduced by imposing the homo- 
condition. 

The method has been generalized to work with multi 
configuration wave functions as in the Multi Configu- 
ration Hartree Fock method. This Multi Configuration 
OEP method provides results very close to those obtained 
by using Multi Configuration Hartree Fock (MCHF). 

An analysis on the convergence on the parameteriza- 
tion of the effective potential is carried out. The results 
obtained are very close to the Hartree-Fock self consis- 
tent energies, eigenvalues and exchange energies improv- 
ing previous optimized effective potential calculations. 
Results for several bound states of the beryllium atom 
by using multi configuration wave functions have been 
reported. 
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